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ABSTRACT 

Poorly understood feedback processes associated with highly-luminous black hole accretion in 
quasars may dramatically affect the properties of their host galaxies. We search for the effect of quasar 
feedback on surrounding gas using Planck maps of the thermal Sunyaev-Zel’dovich effect (tSZ). By 
stacking tSZ Compton-y maps centered on the locations of 26,686 spectroscopic quasars from the 
Sloan Digital Sky Survey, we detect a strong but unresolved tSZ Compton-y signal at >5a signifi¬ 
cance that likely originates from a combination of virialized halo atmosphere gas and quasar feedback 
effects. We show that the feedback contribution to our detected quasar tSZ signal is likely to dominate 
over virialized halo gas by isolating the feedback tSZ component for high- and low-redshift quasars. 
We find that this quasar tSZ signal also scales with black hole mass and bolometric luminosity, all 
consistent with general expectations of quasar feedback. We estimate the mean angularly-integrated 
Compton-y of quasars at z ~ 1.5 to be 3.5x10“® Mpc^, corresponding to mean total thermal energies 
in feedback and virialized halo gas of 1.1 (±0.2) x 10®^ erg, and discuss the implications for quasar 
feedback. If confirmed, the large total thermal feedback energetics we estimate of 5% (±1% statistical 
uncertainty) of the black hole mass will have important implications for the effects of quasar feedback 
on the host galaxy, as well as the surrounding intergalactic medium. 

Subject headings: quasars: general, cosmology: cosmic background radiation 


1. INTRODUCTION 

The coevolutionary view of galaxies and their central 
supermassive black holes (SMBHs), suggested by the 
connection between SMBH masses and the velocity dis- 
persi on of their host galaxy bulg e s (the Mrh-o' relation) 

T-K- 


Ferrarese & Merritt 


2000 


Gebhardt et al. 2000 


Merritt &: h'errarese|200ip , has highlighted the import ant 


role of black hole feedback in t he evolution of massive 
galaxies (for a recent review, see|Fabian|2012|). Feedback 
from the quasar phenomenon is one particular form of ac¬ 
tive galactic nuclei (AGN) feedback, triggered by large- 
scale inflows of gas towards t he central SMBH, possibly 
due to major galaxy mergers dSander s et al.||1988[ |Bah- 
call et al. 1997 Hopkins et al. 20o6 f The parsec-scale 
accretion process onto the blVlBH is radiatively efficient 
and highly luminous, with bolometric luminosities of or¬ 
der erg s ~^, close to the Eddington limit. Radia- 


tive line-driving ([ Murray et al. |1995[ Proga et al. 2000| 
or magnetic forced Rroga|2U03[ Everett|2l)05 ' F'ukumura 


et al.||2014| are belie 
monly observed as 
absorption lines in t 

ved to inevitably drive winds, com- 
alueshifted (Av ~ 10,000 km s“^) 
he ultraviolet (e.g. Weymann et al. 

1991 Gibson et al.|2009 Borguet et ai. 121)131), and in the 

X-rays with near-relativistic blueshitts (e.g 

Reeves et al. 

2009 Chartas et al. 

2009 Tombesi et al. 

)10i2t. 'These 


winds can drive out interstellar gas from the the host 


outflows now frequently obseryed 

Alexander et al. 12010 

Sturm et al.||2011l [Greene et al.||2Ullt [Rupke &: Veilleux 

20ll| 

Greene et al.||2012l [Harrison et al.||2012t [Maiolino 

et al. 

2012 Liu et al. [2013 Harrison et al.[ 2014p. This 
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inflow to their nuclei, th us self-regulating accretion ont o 
their central SMBH (e.g. Silk fc Rees|1998 Fabian|19M ). 

A possible consequence of quasar feedback is the ob- 
served quenching of star formation in the host galaxy 


bulge (Farrah et al. 2012 [Page et al. 2012 Lapi et al. 
20141. This coupling of black hole growth to star tor- 
mation can naturally lead to observed correlations be- 
tween prope rties of the host galaxy and the SMBH mass 


(King 


2003), and can also explain the similar observed 


ity in galaxies (e.g. 

Scannapieco & Oh 

2004 

). Many 

ot these observed trends that may owe t 

0 quasar teed- 


back have been reproduced in semi-analytic models and 
cosmological simulations of galaxy formation that utilize 


yarious sub-grid models to ai 

aproximate quasar feedback 

(e.g. Kauffmann & Haehnelt[ 

2000 iDi Matteo et al. 

2005 

Sprmgel et al. 120051 IBower et 

al.|2006[[Sijacki et al. 

)^007 

Somerville et al.||2008p. Despite these advances, because 


and is often obscured by dust and gas, details of the 
feedback energetics, mechanism, and effects on galaxy 
evolution are still poorly understood. 

Several analytical studies and simulations have sug¬ 
gested that the effects of quasar feedback on the host 
galaxy interstellar medium (ISM) and surrounding cir- 
cumgalactic medium can be probed using the ther¬ 


mal Sunyaev-ZeT dovich effect (tSZ, Sunyaev & Zel- 


dovich|1970l 1972) in the Gosmic Microwave Background 


angular resolution GMB experiments ( 

Voit|1994 Natara- 

[an & Sigurdsson|[1999l [Platania et a 

. 1120021 |Lapi et al. 

2003 

Ghatter.iee & Kosowsky 2007 Scannapieco et al. 

2008 

Ghatterjee et al. 2008p. Since the tSZ effect is a 


electrons, a detection would provide a valuable probe of 
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quasar feedback energetics. However, the quasar feed¬ 
back tSZ signal predicted by these studies is too small 
to be detected for individual quasars. Instead, a cross¬ 
correlation (or stacking) approach between CMB maps 
and large samples of quasars is required with current 
CMB measurements to achieve high signal-to-noise de- 
t ections. This stack i ng app roach is taken here. 


Chatterjee et al. (2010) reported a tentative detec- 
tion of quasar feedback through the tSZ effect by cross- 
correlating photometrically-selected SDSS quasar candi¬ 
dates with WMAP CMB maps; this detected quasar tSZ 
signal was at least an order of magnitude larger than 
that predicted from cosmological simulations tuned to 
match the Mrh-o ' relat ion. However, the origin of the 
Chatterjee et al. (2010) tSZ signal is ambiguous, as it 
may be due to ionized gas in the quasar host galaxies 
heated by quasar feedback and/or virialized halo atmo¬ 
spheres. Indeed, the hot intracluster medium (ICM) in 
massive galaxy clusters, originally heated by virial shocks 
during the cluster formation process, dominates the tSZ 
sky. Aside from the ICM of galaxy clusters, there is 
now evidence from the tSZ effect that virialized halo at¬ 
mospheres also exist in the halos of more moderate-mass 
galaxies. By stacking the tSZ signal from gal axies binned 
in stellar mass using Planck CMB maps, [Planck CoF| 
laboration et al. (2013a hereafter PC13a) report robust 
(>3(7) detection of the tSZ effect in galaxies wi th halo 
masses dow n to log Mhaio = 13.2 Mq (also see Gralla 
et al.|2014D Intriguingly, PC13a find that the ampiitude 
ot this tbZ signal scales as a power-law from clusters 
to more modest-sized galaxies (the Y — Mhaio relation), 
extending from clusters down to galaxies that are nearly 
thr ee orders of magnit ude smaller in halo mass (although 
see 


Greco et al.poid ). 


Separating the quasar feedback and virialized halo 
gas contributions to the quasar tSZ signal can be dif¬ 
ficult without knowledge of the exact quasar host halo 
mass distribution at the redshift of the quasar sam¬ 
ple. Large-scale clustering measurement have established 
that quasars lie in halos with a range of masses with 
mean log Mhaio — 12.4 — 12.9 (e.g. Porciani et al. 

2004 [ Groom et al.||200^ [Myers et al.||^0t>| IShen et al. 


2007 [da Angela et al.[[2008[ [Ross et al.[[2009| White et, 
^ 2012| Font-Ribera et ai.[|2013p. This canonical mean 
quasar halo mass is nearly independent of redshift over 
the redshift range we probe (^ ^ 3), and has only a weak 
dependence on physical parameters such as the b ack hole 
mass, Mrh. or the bolometric lu minosity, Lboi (|Shen et 
^2009 Trainor fc Steidel[[2012[ ). Although the tSZ sig- 
nal from virialized gas in halos of Mq (the fidu¬ 

cial mean quasar halo mass) is expected to be negligi¬ 
ble based on the power-law Y — Mhaio relation reported 
by PG13a, a small (and poorly-constrained) number of 
quasars lying in larger, cluster-mass halos can signifi¬ 
cantly increase the stacked mean quasar tSZ signal, thus 
biasing the inferred quasar feedback energetics. 

Our approach to isolating the feedback tSZ signal is to 
measure the quasar tSZ signal for a subsample of high- 
redshift (z ~ 2 ) quasars. Since cluster-mass halos effec¬ 
tively do not exist at high redshifts, the mean tSZ signal 
for this subsample is not expected to be contaminated 
from the IGM of the small number of cluster-mass ha¬ 
los. To confirm our results, we will also compare our 


high-redshift quasar tSZ signal to that from a separate 
subsample of low-redshift (z ~ 1) quasars. We will show 
that the low-redshift quasar tSZ signal is roughly consis¬ 
tent with that expected from a combination of (1) quasar 
feedback with similar feedback parameters as for high- 
redshift quasars, and (2) a non-negligible tSZ signal from 
virial halo gas that is consistent with estimates based on 
the quasar host halo mass distribution. 

In this paper, we will utilize a tSZ Gompton-y map 
produced using CMB maps from Planck. Aside from 
the greatly improved angular resolution and sensitivity 
in comparison to WMAP^ Planck's large number of fre¬ 
quency channels (six in the High Frequency Instrument) 
is well-suited for tSZ studies. These channels span a 
wide frequency range, including a band centered approx¬ 
imately at the 217 GHz tSZ anisotropy null, as well as 
an 857 GHz band that is dominated by thermal dust 
emission. This enables the spectral separation of signal 
components, and thus high signal-to-noise detections of 
the tSZ effect. 

The outline of this paper is as follows: In Section 2, we 
provide a brief primer on the tSZ effect and quasar blast- 
waves, and then we describe the Planck tSZ Compton- 
y map, dust-correction procedure, and the quasar and 
galaxy samples we utilize. Section 3 describes the tSZ 
map stacking process, the resulting quasar and galaxy 
tSZ signals, and validation of our stacking procedure 
through comparisons to previous Planck results. In Sec¬ 
tion 4, we discuss details of the observed quasar tSZ sig¬ 
nal, the inferred total thermal energetics, dependence on 
quasar physical properties, and implications for quasar 
feedback. We summarize and conclude in Section 5. 
Throughout this paper, we assume a standard ACDM 
cosmology with = 0.272, Ha = 0.728, and Hq = 70.4 
km s“^ Mpc“^, consistent with the WMAP7 results of 


Komatsu et al. (2011). 


2. ANALYSIS METHODS AND DATA 
2.1. The Thermal Sunyaev-ZeTdovich Effect 

In the tSZ effect, inverse Gompton scattering of CMB 
photons by thermal motions of energetic electrons in¬ 
duces a distortion on the blackbody CMB spectrum 
along the line of sight of the form 


AT 

— = y [x coth(a;/2) - 4], 


( 1 ) 


at dimensionless fre quency x = hv/k^T^f (Sunyaev & Zel- 
dovich[[1970[ [1972[), w here Tn is the tempera ture ot the 
GiVIB (for a review, see Carlstrom et al.|2002|. Equation 
I neglects the small relativistic corrections to this spec¬ 
tral distortion, which are only important at extremely 
high temperatures (> 10® K). This spectral distortion is 
parameterized by the dimensionless tSZ Compton-y pa¬ 
rameter integrated along any line of sight I in direction 


n 


y(n) = J UeCTT 


k-eTe 


dL 


( 2 ) 


where and Tg are the number density and temperature 
of the scattering electrons, respectively. 

The total tSZ signal for a source at redshift z is char¬ 
acterized by the angularly integrated Compton-y param- 
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eter 

Yiz) = D\{z) J y{n)dn 


( 3 ) 


integrated over some solid angle D of sky and expressed 
in Mpc^, where Da{z) is the angular diameter distance. 
This integration of the Compton-y over solid angle cap¬ 
tures the integrated Compton-y in a cylinder along the 
line of sight. 

After measuring Y(z) as described by Equation 3, the 
implied total thermal energy in electrons within this ra¬ 
dius is then 

Ee = {z)rfie(? I <7 t . (4) 


For protons and electrons in thermal equilibrium, the 
total energy in ionized gas is simply 

Etot = (1 4 - )Eei ( 5 ) 

Me 

where Me is the mean particle weight per electron, for 
which we will assume the primeval value of 1.17. 

Due to the differences in angular diameter distance, our 
comparison of the mean integrated Compton-y’s both in¬ 
ternally between different redshift bins and to those of 
PC13a (to validate our stacking procedure) is aided by 
a rescaling of the tSZ signal to a common Da = 500 
Mpc. We will aim to show that the detected quasar tSZ 
signal is inconsistent with simply originating from viri- 
alized halo gas, which is the dominant signal measured 
in PC13a. Thus, it is also helpful to divide out the ex¬ 
pected redshift dependence of the integrated Compton-y 
of virialized gas in galaxies, which we can do if we as¬ 
sume self-similar redshift scaling of isothermal virialized 
halo gas in hydrostatic equilibrium. Following PC13a, 
we thus define the rescaled integrated Compton-y: 


Y = 


Y{z) 

D\{z) 


E-^/^{z) 


da{z) y 

500 Mpc ) 


( 6 ) 


expressed in arcmin^, where E{zy = ni„(l + z)'^ + ^a- 
For self-similar scaling and if the virialized halo gas dom¬ 
inates the tSZ, the quantity Y will be constant with red¬ 
shift at fixed halo mass. 

Although there are likely to be correlations between 
quasars and other CMB secondary anisotropies, they do 
not significantly affect our analysis or results. Specifi¬ 
cally, i n the stacking approach we take, the kinetic SZ 
effect (Sunyaev & Zeldovich 1980) will average to zero 
due to the random line of sight velocities of qu asars. Sim¬ 


ilarly, for gravitat ional leasing of the CMB (Blanchard 
& Schneider jl987) by quasar host halos, the sign of the 


fensing anisotropy is dependent on whether hot or cold 
regions of the CMB are lensed, which average to zero in a 
stack. Finally, our estimate of the non-l inear integrated 
Sachs-Wolfe effect (|Rees & Sciama|1968|) shows that it is 
negligible for quasar host halos. 


2.2. Quasar Wind Bubbles 

Quasars are known to drive winds that can be com¬ 
parable energetically to their electromagnetic emissions. 
At the shocked interface between the wind bubble and 
the ambient medium, roughly half of the energ y in the 
outflow is thermalized (Ostriker fc McKee|1988) and can 


hence be visible in a Compton-y map. If the wind bub¬ 
ble is energy conserving, the physics of the bubble would 
depend little on the energy injection process and would 
extend into the circumgalactic medium (and eventually 
the interg alactic medium). The radius of this blastwave 
would be (iKoo & McKee[|1992 1 


Rbw - 0.2 Mpc, (7) 

where Aioo is the mean enclosed pre-shock gas density 
in units of 100 times the cosmological mean, and .Zi 5 = 
(1 -I- 2 :)/ 2 . 5 . Equation 0 assumes a continuous energy 
injection over a time ts with luminosity where 

£'62 is the total injected energy in units of 10®^ erg, and 
tg is the time over which the feedback energy is injected in 
units of 1 0^ yr. Quasar lifetime estim ates favor tg = 0.1 — 
10 (e.g. Martini & Weinberg 2001), and the reference 
energy of 10'^^ erg is 5% of the rest mass enerw of a 
10 ®Mq black hole, a number we will revisit later^ For 
the 10' FWHM resolution of the Planck Compton-y map 
used here (Section 2.3), a distance of {3.7,4.8,5.1} Mpc 
is resolved for 2 : = (0.5,1,1.5}, respectively, distances 
much greater than the expected i?BW while the quasar 
is active. Hence, the Compton-y profile around quasars 
is unlikely to be resolved in this study. 

Equation ([^ is valid if the injected energy is not ra¬ 
diated away and is hence en ergy conserving. [Faucher- 
Giguere & Quataert (2012a) argued that this would be 
the case for quasar winds. In particular, they showed 
that at small radii, where the inverse-Compton cooling 
time is short, the proton-electron equilibration time is 
long enough that the bulk of the energy is carried by the 
protons, making quasar winds energy conserving (at least 
for fast winds with injection velocities >10^ km s“^). In 
order for the wind energy to later produce an observed 
Compton-y signal, the proton bath needs to recouple to 
the electrons well downstream where the cooling time is 
long. We find that this is likely as the electron-proton 
equilibration timescale via Coulomb scattering at i?BW 


IS 


, 3/5 7-9/5,-9/5 A - 3/5 


+ A... c, 

(-ep ~ J^62 ^1.5 '' 


^100 


Myr, 


( 8 ) 


By evaluating tep at the blastwave temperature of 
kbTsw = 3/32MeWpi?^^ and at a post-shock density en¬ 
hancement of four as implied by the Rainkine-Hugoniot 
jump conditions for a strong adiabatic shock. Equation 
8 ignores conduction, which would act to further reduce 
tep. Thus, tep is smaller than the quasar lifetime tg for 
a range of plausible parameters, motivating searches for 
the imprint of quasar feedback through the tSZ effect 
around quasars. 

This picture of energy-driven winds is supported by 
recent investigations of the driving of galactic-scale out¬ 
flows by quasar winds using i dealiz ed and cosmological 
simulations by |Costa et al. (2014), who showed that 
the momentum-tfux from energy-conserving winds are re¬ 
quired to match observed outflows. Furthermore, X-ray 
observations of quasars also do not appear to show signa¬ 
tures of strong inverse-Compton cooling of the shocked 


® If we had instead assumed all of the energy had been injected 
instantaneously a time ig ago (which is the Sedov-Taylor blast- 
wave), the shock radius Ts would be larger by a factor of 2. This 
limit might apply if most fe edback occurred duri ng an earlier time, 
such as an obscured phase ([Hopkins et al.|2008[|. 
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wind (|Bourne & Nayakshin 2013), suggesting that vast 
amounts ot thermal energy lies in the ionized gas sur¬ 
rounding quasars. 


2.3. Planck Compton-y Map 

We utilize the public ly available tSZ Compton-y map 
of Hill & Spergel (2014 hereafter HS14j^ produced using 
15.5-month nominal rriission data irouiPlanck. We only 
briefly summarize its characteristics relevant for our in¬ 
vestigation and refer the reader to HS14 for more details 
on its construction and validation. This Compton-y map 
was created using a modified constrained In ternal Linear 
Combination (ILC, |Remazeilles et al.||2011) method, us¬ 
ing the 100, 143, 217, 353, and 545 GHz channel CMB 
maps from Planck's High Frequency Instrument (HFI). 
In this modified ILC algorithm, the observed emission in 
the channel maps are assumed to be linear combinations 
of blackbody primary CMB emission, tSZ spectral dis¬ 
tortion (parameterized by the tSZ Compton-y in Equa¬ 
tion 1), and noise. Specifically, the CMB thermodynamic 
temperatures Ti{p) in each pixel p of the i channel maps 
are modeled as 


Trip) = divip) + bisip) P Uiip), (9) 

where yip) is the tSZ Compton-y, s(p) is the black- 
body primary CMB, and n(p) is noise. A tSZ Compton- 
y map is generated by computing a set of weights Wi 
that recovers the Compton-y, such that yip) = WiTiip). 
These weights are computed by minimizing the variance 
in the CMB temperature maps and the assumed emis¬ 
sion model, while being constrained to produce zero re¬ 
sponse to the primary CMB emission and unit response 
to a tSZ spectral distortion. The explicit zero response 
to primary CMB emission in this ILC algorithm reduces 
the variance in the resulting stacked maps in compari¬ 
son to stacking of the raw frequency maps, which aids in 
recovery of the low signal-to-noise tSZ effect. We note 
that the assumed emission model does not include dust 
emission, which we will discuss in more detail in Sec¬ 
tion 2.4. The channel maps used in this procedure and 
the final Compton-y map we utilize in our stacking are 
in HEALpi Hg format at a resolution of Nside = 2048, all 
smoothed to a common Gaussian beam with 10' FWHM. 

HS14’s primary departure from the well-known ILC 
algorithm in constructing the final tSZ Compton-y map 
lies in the variance minimization of linear combinations of 
the individual channel maps. Instead of a minimization 
based on the auto-spectra of the channel maps, HS14 uses 
cross-spectra of independent channel maps from Planck's 
single-surveys (‘survey 1’ and ‘survey 2’) in the variance 
minimization. Furthermore, the linear combination vari¬ 
ance is minimized over a restricted multipole range of 
300 < I < 1000, down to scales of approximately the 
beam FWHM. This restriction on the variance minimiza¬ 
tion focuses the minimization on scales of the cosmic in¬ 
frared background and unresolved IR sources rather than 
Galactic dust, which is heavily masked before minimiza¬ 
tion (see below). The final Gompton-y map we use is an 
inverse variance-weighted co-addition of the two single¬ 
survey Gompton-y maps. 


http://www.astro.princeton.edu/ jch/ymapvl/ 
® http://healpix.sourceforge.net/ 



-4.27725e-05 tSZ Compton-y 7.19625e-05 



-4.27725e-05 tSZ Compton-y 7.19625e-05 



-4.64444e-06 tSZ Com'pton-y -2.68542e-07 


Figure 1. Mollweide projection of the tSZ Compton- v map 
(mean-subtracted for visualization) of|Hill &; Spergel|(|2014f, using 
a linear color bar (top), and a high-contrast histogram-normalized 
color bar (middle), as well as the expected dust Compton-y con¬ 
tamination map (bottom). The tSZ Compton-y map used in our 
stacking is corrected for dust contamination by subtracting this 
dust Compton-y contamination map (see Section 2.4). Various 
masks described in Section 2.3, including a dust mask based on 
the Planck 857 GHz channel map and a bright tSZ clusters mask, 
have been applied. 


Before construction of the Gompton-y map, HS14 ap¬ 
plied several masks to each of the channel maps, includ¬ 
ing those to account for the incomplete sky coverage of 
the single-survey maps, all 5a point sources detected in 
any frequency band in both the Low and High Frequency 
Instruments (which removes all resolved radio sources), 
and regions with significant Galactic dust emission. Mas¬ 
sive galaxy clusters are the brightest individual sources 
in the final tSZ map of HSI4, but their presence in our 
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Compton-y map stacking procedure (see Section 3) can 
strongly increase the variance of the background in our 
mean stacked maps, which instead aim to detect tSZ sig¬ 
nal from halos of much lower masses and tSZ Compton- 
y’s. We thus apply an additional bright tSZ cluster mask 
to the final Compton-y map of HS14, based on the Planck 
Early Su nyaev-ZePdovich Detected Cluster Candidates 
Catalog (Planck Collaboration et al. 2011); all pixels 
within the cluster angular extent listed in the catalog 
based on X-ray or Planck observations (where the X-ray 
extent is not available) are masked. 

To minimize dust contamination in their Compton-y 
map, HS14 preemptively masked regions of strong dust 
emission in the channel maps by removing all HEALpix 
pixels with intensity above the 30th percentile in the 
Planck 857 GHz zodi-corrected channel map. However, 
HS14 show that the distribution of the tSZ Compton-y 
in all pixels in their map shows a non-Gaussian tail to¬ 
wards negative Compton-y’s, indicating that some resid¬ 
ual Galactic dust contamination still persists. We thus 
create an additional dust mask using the same method of 
HS14, but instead placing the 857 GHz intensity thresh¬ 
old to remove all pixels above the 20th percentile. In 
our tests, forgoing this additional dust mask and the 
bright tSZ cluster mask increases the variance in the 
background of our stacked Compton-y masks, thus re¬ 
ducing the signal-to-noise of our detections in Section 3. 
The union of our final mask and the mask of HS12 leaves 
a Compton-y map covering 16.4% of the sky (in com¬ 
parison to the 25.2% in the HS14 map), but the map is 
still well-matched to the footprint of the SDSS quasar 
and galaxy sample that we will use. Figure 1 (top panel) 
shows the the resultant masked Compton-y map (mean- 
subtracted for visualization) used in this investigation. 


2.4. Corrections for Dust Contamination 

The primary drawback of using the tSZ Compton-y 
map of HS14 for our analysis is that it does not include a 
dust component in the emission model. When multiplied 
by the ILC weights, dust emission in the five channel 
maps used in the ILC algorithm can produce a non-zero 
response, allowing an additional ‘false’ Compton-y signal 
to leak into the final Compton-y map. Even after impos¬ 
ing our more stringent Galactic dust mask, the Compton- 
y map of HS14 is not completely free of dust contami¬ 
nation. We find that the tSZ Compton-y signal recon¬ 
struction in HS14’s modified ILC algorithm produces a 
negative response to dust emission, primarily due to the 
negative weighting to the high-frequency 545 GHz chan¬ 
nel that is produced by the ILC variance minimization. 
In other words, the ILC weight Wi that is multiplied by 
the thermodynamic temperature in the 545 GHz channel 
map is negative, which can cause the final Gompton-y to 
be under-predicted for strong dust emission in this high- 
frequency band. This can result in systematically under¬ 
predicted or even negative tSZ Compton-y’s in regions 
of strong dust emission, especially for dust emission with 
intensities that increase with frequency. These effects on 
the resulting Compton-y maps are evident in the negative 
Compton-y regions at low Galactic latitudes in Figure 1 
(top panel), and more clearly seen when the color con¬ 
trast is changed from linear to histogram-normalized in 
Figure 1 (middle panel). 

Figure 1 (middle panel) also reveals large-scale fluc¬ 


tuations in the Gompton-y background, which shows 
systematic differences in the mean Compton-y across 
large regions (e.g., the south Galactic cap has system¬ 
atically larger Compton-y than the north). These fluc¬ 
tuations are in part due to the limited multipole range of 
300 < £ < 1000 in the variance minimization of HS14’s 
modified ILC algorithm. This limited multipole range 
in the variance minimization allows dust emission on 
large scales to produce systematic fluctuations in the 
Compton-y map at lower multipoles, which are primarily 
driven by large-scale Galactic dust emission. This issue 
is also evident in the auto-power spectrum of the origi¬ 
nal Compton-y map (see Figure 9 of HS14), where the 
power on scal es i 5 , 300 is systematically higher than 
that found by Plarick Collaboration et ahj (2013b). To 


mitigate the effects of this non-uniform background in 
our analysis, we will subtract the local background in a 
20-30' annulus around each object in our stacking pro¬ 
cedure (see Section 3). Although this background nor¬ 
malization in our stacking minimizes the effects of the 
large-scale Compton-y fluctuations, it does not fully take 
into account the negative Compton-y response in the ILC 
algorithm for intrinsic dust emission from quasars and 
galaxies. This is especially relevant for our investigation 
since strong negative k-corrections of quasar and galaxy 
dust emission in the millimeter/sub-mm spectrum can 
cause the observed dust emission flux densities to be 
nearly redshift-independent. 

To explicitly correct for negative Compton-y response 
in the ILC algorithm from dust emission intrinsic to 
galaxies and quasars, we utilize the 857 GHz Planck 
zodi-corrected nominal channel map (smoothed to the 
common 10' FWHM) as a tracer for the expected dust 
emission in the other channels utilized in HS14’s modified 
ILG algorithm. We assume a fiducial single-temperature 
modified blac kbodv for the dust with spectrum f^, oc 
(e.g. Blain et al.j 2003), where (3 = 2. The 
dust temperature in this fiducial spectrum is fixed to 
Tdust = 20 K and 34 K in the rest-frame for our galaxy 
and quasar stacks, respectively, consistent wit h inferred 
dust tempe ratures in star-forming galax ies (e.g. Clemen s 
et al.|2013 ) and quasars at z ~ 1.5 (e.g. Dai et aC |2012| . 
4'he redshifted dust temperatures 7dust/(l -I- z) at the 
mean redshifts of our galaxy and quasar samples are con¬ 
veniently all ~13.3 K, and so we adopt this single red- 
shifted dust temperature for our fiducial dust spectrum. 

Using the above fiducial dust spectrum with redshifted 
Udust/(1 + z) = 13.3 K, we extrapolate the brightness 
temperature of dust emission in each of the other HFI 
channels based on the normalization provided by the 857 
GHz channel map and the effective frequencies of the 
HFI channels provided in the Planck Explanatory Sup- 
plemenl[^ We then produce a map of the expected neg¬ 
ative Gompton-y response due to dust emission by mul¬ 
tiplying the expected dust emission in each channel (in 
thermodynamic temperature units) by the ILC weights 
of HS14 for the public 30% threshold map we use (see 
their Table 2). The resulting dust Compton-y map is also 
shown in Figure 1 (bottom panel), and shows systemat¬ 
ically negative Compton-y’s due to dust contamination. 
A stacking test of our dust map on our quasar sample 
shows an excess centered on the quasars that corresponds 


® http://www.sciops.esa.int/wikiSI/planckpla 
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Figure 2. Redshift distributions for our quasar and Locally 
Brightest Galaxy samples. The median redshifts of the quasar, 
lower- and higher-mass LBG samples are 2 = 1.5, 0.49, and 0.54, 
respectively. 

to integrated dust luminosities of order erg s“^. To 
correct HS14’s tSZ Compton-y map for this dust emission 
intrinsic to galaxies and quasars, we subtract this dust 
Compton-y map (thus increasing the Compton-y’s), and 
the final dust-corrected tSZ Compton-y map is used in 
our quasar and galaxy stacking procedure (Section 3). 

We empirically tested the sensitivity of our results to 
the assumed dust spectrum by also producing stacked 
Compton-y maps for our quasar and galaxy samples us¬ 
ing various spectral indices /3 from 1.0 to 2.0 (which 
spans the range exp ected for the rmal emission in the¬ 
oretical dust models; |Draine||2011|), as well as redshifted 
dust temperatures Tdust/(1 + z) irom 5 K to the high- 
temperature Rayleigh-Jeans tail. The resulting differ¬ 
ence in the Compton-y’s is ^10% and are stable around 
our fiducial dust spectrum, showing that our results are 
robust to the assumed dust spectrum. 

2 . 5 . SDSS Spectroscopic Quasar Sample 

We use a sample of spectroscopically identi fied quasars 
from the Sloan Digital Sky Survey (SDSS, [ York et ^ 
20001) Data Release 7 (DR7) Quasar Catalog ([Schneider 
et al.|2007[) We also utilize various additional pr operties 
of t hese D R 7 quasars calculated and complied in Shen et 
al. (20071, including virial black hole masses Mbh, bolo- 
metric luminosities Tboh and radio flux densities. This 
DR7 quasar catalog contai ns 105,783 quasars, s elected 
primarily by optical color (Richards et al. 2002). The 
spectroscopic footprint covers approximately 10,000 deg^ 
over the north and south Galactic caps, largely overlap¬ 
ping the Planck tSZ Compton-y map. From this DR7 
catalog, we select 39,117 quasars which fall in HEALpix 
pixels not masked in the various tSZ Compton-y map 
masks described in Section 2.3. 

The tSZ effect at millimeter wavelengths from quasars 
may be susceptible to contamination from intrinsic ra¬ 
dio emission, particularly since approximatel y 10% of 
spectroscopic quasars in SPSS are radio-loud (Ivezic et 


al. 2002 Balokovic et al. 2012). Although the Planck 


point sources mask imposed on the Compton-y map re¬ 
moves all resolved radio sources detected in any Planck 
frequency channel, some unresolved radio emission may 
still be present. We thus remove all quasars in our sample 
with observed radio flux density at rest-fram e 6 cm above 
2 mJy observed by the F IRST radio survey (Becker et al. 
1995 White et al.||1997), as well as a small number out- 


coadd tSZ Compton-y map 
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Figure 3. Mean tSZ map from our stack of 26,686 quasars. A 
clear tSZ signal centered on the quasars is evident. The coherence 
length of structures in this stacked map is set by the 10' FWHM 
Gaussian beam, and the quasar tSZ signal is unresolved. 

side the FIRST survey footprint for which the radio flux 
densities are unconstrained (in total ~10% of our sam¬ 
ple) . A separate stacking test of our quasar sample that 
includes the radio-bright quasars in our sample produces 
only slightly lower mean Compton-y’s. Furthermore, we 
find that even if all of our radio-quiet quasar sample had 
fluxes just below the FIRST point source threshold, they 
would contribute just 20% to the Y we detect (ass uming 
a mean spe ctral index in specific intensity of — 0.7; [Ivezic| 
et al.|2002 ). Thus, we conclude that radio emission does 
not significantly contaminate our stacks. 

2.6. SDSS Spectroscopic Locally Brightest Galaxy 
Sample 

To verify that the mean quasar tSZ signal measured 
from our Compton-y map stacking procedure is robust, 
we also stack samples of galaxies and compare their tSZ 
signal to those measured by PC13a. Specifically, we will 
use SDSS spectroscopic galaxies and construct samples 
of Locally Brightest Galaxies (LBGs), following the pro¬ 
cedure of PC13a. These LBGs are the most luminous in 
their local vicinity, and are thus likely the main galaxy in 
their halos. The unobservable halo mass of each LBG in 
our sample can be estimated using estimates of their stel¬ 
lar mass, along with the stellar-to-halo mass (SHM) re¬ 
lation calculated by PC13a. This SHM relation is based 
on an updated ver sion of the semi-an alytic galaxy forma¬ 
tion simulation of Guo et al. (2011), which PG13a used 
to select a simulated LBG sample, and then calculated a 
SHM relation based on their stellar and halo masses in 
the simulation. 


We use the SDSS Data Release 10 (DRIO, Ahn et al. 
2014) spectroscopic galaxy sample, which includes stel¬ 
lar iiiass estimates from the Portsmouth Group based 
on their SDSS ugriz photometry. These galaxy stellar 
ma sses a re estimated using the method of|Maraston et 


al. 


et al. 


an d the stellar populations model ot Maraston 
(|2009l). We select spectroscopic galaxies with ob- 
served tab < 21 model mag from the photometry, and 
spectroscopic redshifts z > 0.03. To remove satellite 
galaxies, we select samples of LBGs following the method 



















































































Detection of Quasar Feedback 


7 


of PC13a, in which the SDSS photometric galaxy cata¬ 
log is used to remove spectroscopic galaxies with brighter 
nearby photometric companion galaxies. Specifically, we 
remove any spectroscopic galaxy that has a nearby pho¬ 
tometric galaxy with brighter r model mag, and is within 
a projected distance of 1.5 h~^ Mpc and redshift sepa¬ 
ration <2000 km s“^. The redshifts of the photometric 
galaxies used in this procedure are photometric redshifts, 
es timated using the kd-t ree nearest neighbor fit described 


Csabai et al. (20031 and provided as part of SDSS 


DKIO. This selection of LBGs allows us to directly use 
the SHM relation of PC13a to estimate their halo masses, 
and it was shown in PC13a to be robust against changes 
in the projected distance parameters of 1.5 h~^ Mpc and 
2000 km s“^, with low contamination (< a few percent). 

We first construct a lower-mass LBG sample with halo 
masses similar to the mean quasar Mhaio = 10 ^^'^ Mq; 
these are approximately the lowest mass galaxies for 
which tSZ signal is still marginally detected by PC13a. 
Specifically, we select a sample of 70,889 LBGs with 
stellar masses 10 ^° ® < M* < 10 ^^'^ Mq that are not 
in masked pixels of our tSZ Compton-y map. LBGs 
with these stellar masses lie in halos with mass approx¬ 
imately in the range 10^^-® < M^aio < 10^^ ® Mq. We 
note that halo masses provided by the SHM relation of 
PG13a are presented in units of M 500 , while the observed 
quasar halo masses estimated from clustering measure¬ 
ments we compare to are ge nerally in units of M -im- For 
Navarro-Frenk-White halos ( Navarro et al. lip'l l of mass 
-^200 = Mq (similar to quasars), and a concentra¬ 

tion p arameter consistent with that found by Neto et al.| 
(2007), M 500 is 10^^ ® Mq. Since M 500 and M 200 are not 
significantly different at the halo mass scales we probe, 
we treat them as equivalent quantities and refer to them 
as M\iqIq . 

We also construct a sample of higher-mass LBGs that 
have stronger tSZ signal, for comparison to PG13a. We 
select 81,766 galaxies with stellar mass < M^, < 

10^1® Mq (corresponding to halo masses of 10^® ® < 
A7haio < 10^®-® Mq). The details of these two LBG sam¬ 
ples are also summarized in Table 1, and their redshift 
distributions are also shown in Figure 2, with median 
redshift z = 0.49 and 0.54 for the lower- and higher-mass 
LBG samples, respectively. 


3. STACKING COMPTON-Y MAPS OF THE QUASAR AND 
GALAXY SAMPLES 

To detect tSZ signal from our quasar and Locally 
Brightest Galaxies (LBGs) samples, we will stack regions 
of the Gompton-y map described in Section 2.3, centered 
on each of the different samples. We will first do this in 
angular space for an initial detection, and then in phys¬ 
ical space for the actual measurement of the mean inte¬ 
grated Compton-y’s. While PC13a also showed stacks of 
the tSZ (constructed using an ILG algorithm) centered 
on their LBG galaxy samples, their actual tSZ measure¬ 
ments were made using a multi-frequency matched filter 
(MMF) algorithm. The MMF algorithm used by PG13a 
assumes a template tSZ profile for the sources, perhaps 
better-suited for the tSZ signal from galaxies where there 
is some understanding of the pressure prohle. However, 
the pressure profiles of ionized gas in quasar host halos 
due to feedback is unconstrained. Therefore, it makes 
sense to stack around the positions of the quasars, and 


apply an aperture photometry approach to empirically 
measure the integrated Gompton-y’s without assuming 
a radial template. We will verify that this approach on 
our LBG samples yields a signal that is consistent with 
those reported in PG13a. 

3.1. Angular-Space Stacking 

For each quasar in our sample, we cut out a 1° by 
1° square region from the masked tSZ Compton-y map, 
projected in a gnomonic projection and centered on the 
quasar. Each of these regions is comprised of 60x60 pix¬ 
els at resolution I'/pixel, and their rotational orientation 
are chosen at random. We will not include any quasars 
in our stack for which more than 10 % of pixels within 
30' of the central quasar are masked. To mitigate the ef¬ 
fects of large-scale fluctuations in the background of the 
tSZ Compton-y map noted in Section 2.4, we subtract 
the mean Compton-y in a background annulus of 20-30' 
around each quasar. We produce an equal-weight mean 
co-added map of the remaining 26,686 quasars in our 
sample, which is shown in Figure 3 and has an implicit 
10' FWHM Gaussian beam from the original Compton-y 
map. A clear excess is present in this stacked map peak¬ 
ing at the location of the quasars. The above process is 
repeated to also create separate co-additions of our LBG 
samples. 

Figure 4 shows the radial profile of the tSZ Compton- 
y signal from the stacked image of our quasar sample, 
beginning from the center of the image in Figure 3, using 
bins of width 1.25'. Figure 4 also shows the Compton-y 
radial profile of a stacked quasar map without our dust 
corrections for comparison; as expected, correcting for 
the negative ILG Compton-y response to intrinsic dust 
emission increases the tSZ signal. We perform a null 
test of our detected quasar tSZ signal by stacking on 
the locations of 2 x 10 ® randomly-distributed points on 
the Compton-y map for which <10% of the pixels in the 
stacked region are masked. The Compton-y profile of this 
stacked map is also shown in Figure 4, and is consistent 
with a non-detection. A separate null test (not shown) 
in which the position of each quasar in our sample is 
randomly offset by 1° in Galactic latitude before stacking 
also shows a similar non-detection. 

We produce separate stacked tSZ images of our two 
LBG samples (summarized in Table 1), and their mean 
Compton-y radial profiles are shown in Figure 5, along 
with the dust-corrected quasar radial profile from Fig¬ 
ure 4 for comparison. The lower-mass LBG sample does 
not show a significant tSZ signal, while our higher-mass 
LBG sample shows a clear detection; these results are 
qualitatively consistent with those of PC13s for LBGs in 
these mass ranges, and we make a quantitative compar¬ 
ison of their integrated Compton-y’s in Section 4.1. We 
note that the samples of quasars and LBGs used in the 
stacked Compton-y maps in Figures 4 and 5 are smaller 
than the original samples defined in Sections 2.5 and 2.6 
due to the requirement that < 10 % of pixels in the stacked 
regions are unmasked. The la uncertainties on the radial 
profiles in Figures 4 and 5 are estimated using a standard 
bootstrap resampling of the individual regions compris¬ 
ing the stacked map. Since the number of pixels in each 
azimuthal bin increases with radius, the uncertainties in 
the Compton-y radial profiles are naturally smaller at 
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Table 1 

integrated Compton-y values for our samples of quasars and locally brightest galaxies 


Sample 

Halo Mass 
[Mq] 

Median 

Redshift 

•ya, 

[10~^ arcmin^] 

yb 

[10-® Mpc2] 

pc 

■Alot 

[10®^ erg] 

SDSS DR7 Quasars 

median logigM^aio — 12.7 

1.50 

96 ± 19 

3.5 ± 0.7 

11 ± 2.2 

SDSS DR7 Quasars, higher -2 

median logigM^aio — 12.7 

1.96 

115 ± 19 

4.8 ± 0.8 

15 ± 2.8 

SDSS DR7 Quasars, lower -2 

median logioAfhalo - 12.7 

0.96 

74 ± 30 

2.2 ± 0.9 

7.1 ± 3.1 

LEG, 10.9 < logioM* < 11.1 

12.5 < log^gMhaio < 12.8 Mq 

0.49 

4 ± 20 

0.1 ± 0.5 

0.3 ± 1.6 

LBG, 11.4 < logj^pM* < 11.6 

13.6 < logioMhalo < 13.8 Mq 

0.54 

56 ± 16 

1.4 ± 0.4 

4.5 ± 1.3 

virialized quasar halo gas, 5-param HOD 


0.5 

42 

1.0<^ 

3.2'i 

virialized quasar halo gas, 6-param HOD 


0.5 

27 

0.7^ 

2.3*^ 


^ Redshift-independent integrated Compton-y, Equation 6 
^ Integrated Compton-y, Equation 3 
^ Total thermal energy in ionized gas, Equation 5 

^ These Y and Etot values for virialized halo gas in quasar host halos at 2 : = 0.5 have been rescaled to the 2 : = 0.96 median redshift 
of our low -2 quasar sample for direct comparison. See Section 4.1 for details. 




arcmin from center 


arcmin from center 


Figure 4. Radial Compton-y profile of the mean stacked map 
centered on the quasar sample with correction for dust contami¬ 
nation (red line), and without dust corrections (blue line), as well 
as a null test (green line). Dust contamination in the ILC al¬ 
gorithm produces a negative Compton-y response that we find is 
robust to the dust spectral template, and thus the quasar tSZ sig¬ 
nal is stronger after dust corrections. The Icr uncertainties shown 
(shaded regions) are estimated using a bootstrap resampling of the 
quasar-centered cutouts comprising each stacked map. 

large radii. 

3.2. Physical-Space Stacking 

Due to the different angular diameter distances of our 
Locally Brightest Galaxies (LBGs) and quasar samples 
(as well between objects in each sample), their mean in¬ 
tegrated Compton-y (Equation 3) and thermal energy 
in electrons (Equation 4) should be measured and com¬ 
pared in physical space (although we note that the angu¬ 
lar diameter distance only varies by approximately 50% 
between z = 0.5 and z = 2, making the angular and 
physical space stacks qualitatively similar). We produce 
a separate stacked image of the quasar sample using re¬ 
gions of 30 h~^ Mpc by 30 h~^ Mpc in physical distance 
at the redshift of each quasar. At the median quasar 
redshift of z = 1.5, this is approximately the same re¬ 
gion in physical space as the 1° by 1° regions used in 
the angular-space stacked image of Figure 3. Intention¬ 
ally scaling the size of the regions in this manner causes 
the resulting physical-space stacked image (not shown) 
to look similar to that in angular-space in Figure 3, and 


Figure 5. Radial Compton-y profile of mean stacked maps (all 
dust corrected) centered on the quasar sample (red line), the lower- 
mass LEG galaxy sample (green line), and higher-mass LEG sam¬ 
ple (blue line). The mean integrated Compton-y of our LEG sam¬ 
ples are consistent with the results of PC13a, indicative of hot 
halo atmospheres in these galaxies. However, the Compton-Y of 
the quasar sample is significantly larger than that expected from 
halos of similar masses (compare green and red curves), suggesting 
additional heating from quasar feedback (see Section 4.1). 

allows US to use a background annulus of 10-15 h~^ Mpc 
around each quasar for background subtraction (similar 
to the 20-30' in the angular-space stacks). Similarly, for 
the physical space stacks of our LBG samples, we use 
regions of 20 h~^ Mpc by 20 h~^ Mpc at the redshift of 
each LBG (also approximately a 1° by 1° region at the 
LBG median redshift of 0.52), and background annulus 
of 6.67-10 h~^ Mpc. The similarity in size of the stacked 
regions in physical and angular spaces also causes the 
Compton-y radial profiles of our quasar and LBG sam¬ 
ples in physical space to be similar to those in angular 
space (Figures 4 and 5). 

3.3. Integrated Compton-y Measurements: 

Comparison to PC 13a 

We calculate the mean Y (Equation 3) and Y (Equa¬ 
tion 5) of our quasar and LBGs from their median 
Compton-y radial profiles of our physical space stacks. 
For the integrated Compton-y, we take a simple aperture 
photometry approach to the radial profiles of the previ- 
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ously background-subtracted Compton-y stacked maps 
by integrating their Compton-y profiles within a central 
circular source region. For the quasar and LEG samples, 
this source aperture is chosen to be at radius 5 h~^ Mpc 
and 3.75 h~^ Mpc, respectively. These aperture radii in 
physical space correspond to approximately 10 ' for the 
quasar and LEG samples at their median redshifts in the 
angular-space Compton-y profiles in Figure 4, and twice 
the beam FWHM in the Compton-y map. 

The Y (Equation 5), rescaled angularly and in redshift, 
is estimated based on the resulting Y values, by fixing 
the angular diameter distance to 500 Mpc and assum¬ 
ing the median redshifts of z = 0.52 and 1.5 for our LEG 
and quasar samples, respectively. The measured mean Y 
and Y values of our quasar and LEG samples are listed 
in Table 1. The Icr uncertainties on these integrated 
Compton-y’s are also based on bootstrapping of the in¬ 
dividual regions that comprise their stacked images. 

We verify our tSZ measurements from our stacking 
procedure by comparing our Y measurements for LEGs 
to those of PC13a. Comparison of Y to previous work 
should account for any differences in the solid angle over 
which the Compton-y is integrated, as well as whether 
the Compton-y measurements capture the integrated 
Compton-y in a cylinder along the line of sight or over 
an enclosed sphere assuming a pressure profile. Our mea¬ 
surement of Y for galaxies using Equation 3 is performed 
by integrating the line of sight Compton-y over a solid 
angle corresponding to the 3.75 h~^ Mpc radius of the 
source aperture; this captures the integrated Compton-y 
in a cylinder along the line of sight. In contrast, PC13a’s 
Y measurement using the MMP approach captures the 
integrated Compton-y inside a spherical volume with 
mean density of bQQpc{z), where Pc{z) is the critical den¬ 
sity of the Universe. The radius of this volume (i? 5 oo) 
for galaxies of each halo mass is estimated using the uni¬ 
versal pressure profile of Arnaud et al. (2010), and we 
denote this resulting integrated Compton-y reported by 

PC13aasfX- 

To account for the differences in source aperture sizes 
in our comparison of of LEGs measured by PC 13a 
to the Y we measure, we calculate the conversion fac¬ 
tor Y /Eflioo fo’' galaxy samples. The 3.75 h~^ Mpc 
radius of the circular aperture we use to integrate the 
Compton-y over solid angle for LEGs is approximately 
15.8i?5oo and 6 . 8 i? 5 oo for our lower- and higher-mass 
LEG samples, res pectively. Ey integrat ing the universal 
pressure profile of|Arnaud et al.|(|2010| out to 15.8i?5oo 
and 6 . 8 i? 5 oo in a cylinder along the line of sight (to find 
Y), and dividing by the integrated pressure profile in a 
spherical volume out to R 500 (to find 

Y /Uflfoo = 1. 52. This value of Y /T^foo (consistent with 
that found by Greco et al.|2014 ) is the same for both our 
lower- and higher-mass LEG samples despite the differ¬ 
ent aperture sizes we use of 15.8i?5oo and 6 . 8 i? 5 oo, respec¬ 
tively, because the gas pressure in the universal pressure 
profile is negligible beyond 5 i? 5 oo- 

Por our lower-mass LEG sample (with 10^° ® < M* < 
Mq), our measured Y = 4(±20) x 10“® arcmin^ 
is consistent with the Y = (1.9 ± 1.2 to 2.3 ± 0.9) x 10“® 
arcmin^ detected at low significance by PC13a (see their 


Table 1, multipli ed by Y /Y^^jL = 1 .52), as well as the 
null detection by Greco et al. ( 2014[ ). The two values of 
Y measured by PC13a that we are comparing to is due to 
their finer binning of LEG samples in stellar mass, where 
each of our M* bins encompasses two bins in PC13a. Eor 
our higher-mass LEG sample (with 10^^'^ < M* < 10^^ ® 
Mq), the F = 56 (±16) x 10“® arcmin^ we measure is in 
excellent agreement with the strong Y = (44±5.8 to 91± 
10) X 10“® arcmin^ measured by PC13a. This recovery 
of the tSZ signal from LEGs confirms that our stacking 
method recovers similar tSZ signals as PClSa’s MMF 
approach. 

4. A FEEDBACK ORIGIN FOR THE QUASAR TSZ SIGNAL 

4.1. Separating Quasar Feedback 
and Virialized Halo Gas 

The quasar tSZ signal we detected in our Compton-y 
map stacks is likely to be a combination of both feedback 
effects and virialized halo gas in quasar host galaxies. It 
may seem that the feedback effects can be isolated by 
simply taking the total quasar tSZ signal and subtract¬ 
ing off the small tSZ signal for Locally Erightest Galax¬ 
ies (LEGs) with halo masses similar to the 10^^'^ Mq 
mean halo masses of quasars. However, a more careful 
separation of these two components of the total quasar 
tSZ signal should take into account (1) the full distri¬ 
bution of quasar host halo masses, and ( 2 ) the different 
median redshifts between our quasar and LEG samples. 
Specifically, an estimate of the virialized halo component 
in the quasar tSZ signal based on the 10^^'^ Mq halo 
mass LEGs is likely to be biased due to a small number 
of quasars lying in cluster-mass halos. The tSZ signal 
of these handful of quasars lying in clusters are domi¬ 
nated by the ICM rather than feedback, and are orders 
of magnitude larger than the tSZ signal of more typical 
quasars due to the power-law relation between the tSZ 
integrated Compton-y and halo masses (the Y — Mhaio 
relation). These quasars can thus significantly increase 
the mean quasar tSZ signal in our stacks, and accounting 
for this effect requires knowledge of the quasar host halo 
mass distribution. Furthermore, the number of cluster- 
mass quasar host halos is redshift-dependent, and so the 
quasar tSZ signal should be compared to estimates of the 
virialized halo gas contribution at similar redshifts. 

An approach we take for isolating the quasar feedback 
tSZ signal is to create higher- and lower-redshift quasar 
subsamples, with which we produce separate stacks to 
measure the integrated Compton-y. We argue that the 
high-redshift subsample will not contain quasars lying in 
massive clusters, and thus its tSZ signal will be purely 
due to feedback with little to no contamination from viri¬ 
alized halo gas. In addition, the lower-redshift subsample 
will be used as a verification of the feedback energetics 
derived from the high-redshift subsample. Specifically, 
we will estimate the virialized halo gas component in 
our low-redshift quasar tSZ signal, based on the quasar 
halo mass distribution at z cz 0.5 inferred fr om quasar- 
galaxy cross-correlation by Shen et al. (2013). Since the 


quasar feedback tSZ signal is a measurement of the total 
thermal feedback energy in ionized gas, we will parame¬ 
terize this quasar feedback energy Efb (calculated using 
Equation 5) as a function of the rest-mass energy of the 
accreted material (i.e. Mbh), such that Efb = pcMbhc^- 
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quasar log(JV4,J[Mg]) 

Figure 6. Contribution to the mean tSZ Y (Equation 6) from 
virialized gas in quasar host halos as a function of halo mass, based 
on the halo mass probability distribution from a 5-parameter (red 
solid line) and 6-parameter (blue dashed line) halo occupation dis¬ 
tribution model. 


The key feedback parameter is the observationally de¬ 
generate product tye, where rj is the radiative efficiency 
of accretion, and e is the feedback coupling efficiency 
of radiation to the thermal energy in ionized gas. The 
spectroscopically-derived estimates of Mbh we use for 
each quasar in our sample are single-epoch Mbh esti¬ 
mates based on the broad emission lines, and have pos- 


constrained (iDenney et al. 2009 

Rafiee & Hall 2011 

Den- 

ney et af.||2013p. 


ples, we split our quasars at z ~ 1.5, approximately 
the median redshift of the main quasar sample. The 
high- and low-redshift subsamples have median redshifts 
of 1.96 and 0.96, and median Mbh of 10®'^ and 10® ® 
Mq, respectively. A physical-space Compton-y stack of 
the high-redshift quasar subsample following the proce¬ 
dure in Section 3.2 results in an integrated Compton- 
y of 4.8 (±0.8) X 10“® Mpc^, implying a total thermal 
energy of 15 (±2.8) x 10®^ erg and a quasar feedback 
efficiency from accretion of rje = 0.05 (±0.01 statisti¬ 
cal uncertainty), or approximately 5% of the black hole 
mass. This feedback energy (and the implied feedback 
efficiency) for high-redshift quasars is not significantly 
affected by virialized halo gas, since the mean tSZ sig¬ 
nal for virialized gas in quasar halos at the mean quasar 
halo mass of 10^^-^ Mp , is negligible (as found by PC 13a 
and IGreco et al 




__ and confirmed in our lower-mass 

LEG stack in Section 3.3). More importantly, the num¬ 
ber of quasars lying in cluster-mass halos is also effec¬ 
tively zero at z ~ 1.96 since they will not have yet formed 
at these redshifts. We note that our quoted uncertain¬ 
ties on Tje are statistical only, and do not include the 
poorly-understood systematics in the single-epoch Mbh 
estimates we use. 

We further verify this measurement of the quasar feed¬ 
back efficiency from our high-redshift quasar sample, by 
comparing to an analysis of our low-redshift quasar sam¬ 
ple. At low-redshifts, the non-negligible tSZ signal from 
virialized halo gas must be subtracted to isolate the 
quasar feedback energy. To estimate this contribution. 


we utilize the quasar number density as a functio n of halo 
mass at z ~ 0.5 calculated by Shen et al. (2013|). This is 
inferred from halo occupation distributioii (HUD) model¬ 
ing of the quasar-galaxy two-point cross-correlation and 
when normalized, gives the probability distribution of 
quasar host halo masses. Different quasar HOD models 
are often k nown to be degenerate with clustering mea- 
surements ( Richardson et al.||2012 Kayo & Oguri 2012 


Ghatterjee et al. 2013|), and the uncertainties on the best- 


tit occupation functions are often large. Nevertheless, 
they provide a useful check of our high-redshift feedback 
results, and so we utilize the quasar M^aio distribution 


inferred in Shen et al. (2013) using both a 5-parameter 
and 6-paraineter HUD model (see their Figures 14 e and 
f). We multiply these two quasar Mhaio probability dis¬ 
tributions by the redshift-independent power-law rela¬ 
tion between and Mhaio found in PG13a (see their 

Equations 1 and 2), which results in the contribution to 
the mean tSZ signal from virialized gas in quasar host 
halos as a function of Mhaio at z ~ 0.5, shown in Figure 
6. Since we are comparing results from PC13a to our 
measurements, the integrated Gompton-y measurements 
of PC13a have been multiplied by the extra factor of 


Y /K^ 5 oo = 1.52 to account for differences in the measure¬ 
ment techniques, as discussed in Section 3.3. Although 
cluster-sized quasar host halos are rare, Figure 6 shows 
that halos of ^10^^ ® Mq are the primary contributors 
to the mean virialized halo gas tSZ signal in quasar host 
halos at z ~ 0.5. Since the virialized halo gas (i.e. ICM) 
in these cluster-mass halos are generally well-modeled 
as isothermal gas in hydrostatic equilibrium, our direct 
comparison of the virialized gas contribution in LBGs to 
quasars with similar halo masses should not be affected 
by additional sources of tSZ signal (e.g. stellar feedback 
and non-thermal pressure support). 

Integrating the Y contributions in Figure 6 over all 
halo masses gives the mean Y from virialized gas in 
quasar host halos at z ~ 0.5; for the 5- and 6-parameter 
models, we find Y = 4.2 x 10“® arcmin^ and Y = 
2.7 X 10“® arcmin^, respectively. Since we will be com¬ 
paring this tSZ signal at z ~ 0.5 to our low-redshift 
quasars at z ~ 0.96, the redshift evolution of this viri¬ 
alized halo gas tSZ signal and the difference in angu¬ 
lar diameter distances should be taken into account. 
We thus rescale these mean virialized halo gas Y val¬ 
ues to F(z = 0.96) to find Y{z = 0.96) = 1.0 x 10“® 
Mpc^ and T(z = 0.96) = 0.7 x 10“® Mpc^ for the 
5- and 6-parameter HOD models, also summarized in 
Table 1. These expected mean quasar halo gas tSZ 
Y{z = 0.96) are smaller than our observed low-redshift 
quasar Y = 2.2 (±0.9) x 10“® Mpc^ by approximately 
2 to 3 times for the 5- and 6-parameter HOD models, 
respectively. This suggests the additional presence of 
feedback effects in our low-redshift quasar sample, which 
likely dominates our observed quasar tSZ signal over the 
contribution from virialized halo gas. 

To summarize our method above, we have taken the 
quasar halo mass distribution at z = 0.5, multiplied by 
the Y- Mhaio relation for LBGs, and integrated to find 
the virialized halo gas contribution in our low-redshift 
quasar tSZ Y signal. Since the low-redshift quasars we 
wish to compare to have median redshift of z = 0.96, 






































Detection of Quasar Feedback 


11 


we rescaled the estimated virialized halo gas Y (which is 
redshift-independent) to the redshift-dependent Y{z = 
0.96). However, we note that despite this rescaling of 
Y to Y(z = 0.96), Y(z = 0.96) is still not a completely 
accurate estimate of the virialized halo gas tSZ signal 
for quasars at z = 0.96, because the original quasar halo 
mass distribution we use is instead measured for z = 0.5. 
This will lead us to slightly overestimate the virialized 
halo gas tSZ Y{z = 0.96) by overestimating the num¬ 
ber of cluster-sized quasar host halos at z ~ 0.96, and 
thus underestimate the feedback signal. Limiting our 
low-redshift quasar sample to only quasars at z ~ 0.5 re¬ 
sults in a subsample that is too small to produce a high 
signal to noise stacked tSZ detection. Nevertheless, as 
a useful additional check of the high-redshift feedback 
parameters, we isolate the feedback tSZ signal for our 
low-redshift sample at z ~ 0.96 by subtracting our best 
estimate of the mean virialized halo gas component from 
the 5- and 6-parameter HOD models. We find that the 
isolated low-redshift quasar feedback has mean energies 
of 4.4 (±3.1) X 10®^ erg, giving a quasar feedback effi¬ 
ciency from accretion of rje ~ 0.038 (± 0.027 statistical 
uncertainty). Considering that this is likely to be a small 
underestimate of rye, we conclude that this is crudely con¬ 
sistent with the rye = 0.05 feedback efficiency we measure 
at high redshifts. 


4.2. Quasar tSZ Trends with Mbh, Lhoi, and z 

In Section 4.1, we showed that the effects of quasar 
feedback likely dominate our observed quasar tSZ signal 
at both low and high redshifts. As a test of this quasar 
feedback interpretation, we investigate the dependence 
of our detected quasar tSZ signal on physical properties 
of the quasars, including Mbh, l^boi, and z. Because the 
cooling time of the shocked wind is long, the feedback 
energy Efb we measure is likely a total thermal feedback 
energy released into the surrounding gas over the quasar 
lifetime, and so a positive scaling of t he tSZ signal with 
Mrh and Lhni is generically expected (Scannapieco fc O h 
2004| iScamiapieco et al.||2008| |Chatterjee et al.||2UU8p. 
We divide our quasar sample in two bins of Mbh and 
Lbob and compare their tSZ Compton-y radial profiles in 
separate angular-scale stacked maps of each subsample. 
The spectroscopically-der ived estimates of M rh and Lboi 
we use are provided from Shen et al. (20111. 

We divide our sample into two bins m Mrh (7 < 
logioMRH < 9 Mq, and 9 < log;iQMBH < 10.5 Mq), 
as well as in Lboi (44.5 < log^ol^boi < 46.3 erg s“^, 
and 46.3 < log^oAboi < 48 erg s“^). These bins are 
split at approximately the median Mrh and Lboi, with 
median logj^QAfRR = {8.7,9.3} and median logj^pTboi = 
{45.9, 46.7} in the bins, respectively; this results in simi¬ 
lar signal-to-noise for each of their separate stacked maps 
(all dust-corrected). Figure 7 shows the Compton-y ra¬ 
dial profiles of quasars with different Mrh (top panel) 
and Lboi (middle panel); a clear trend of increasing 
Compton-y signal with both Mrh and Lboi is evident, 
as expected if the observed quasar tSZ signal is domi¬ 
nated by quasar feedback. This feedback interpretation 
is also supported by the positive scaling with Lboi, al¬ 
though time-variability of the quasar continuum emission 
for each quasar will increase scatter in this trend. 

Figure 7 (bottom panel) also shows tSZ Compton- 



arcmin from center 



arcmin from center 


Figure 7. Radial tSZ Compton-y profiles of mean stacked maps 
of quasars (all dust-corrected), separating the stacks into a low and 
high black hole mass sample (top panel), a low and high bolomet- 
ric luminosity sample (middle panel), and a low and high redshift 
sample (bottom panel). A clear trend of higher tSZ Compton-y 
with increasing black hole mass and bolometric luminosities is evi¬ 
dent. The quasar tSZ signal is even somewhat stronger in the high 
redshift sample. In our feedback picture, this trend occurs because 
more massive black holes exist in our sample at high redshifts, and 
it also indicates that the signal cannot be dominated by virialized 
gas of cluster-mass halos. 
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y radial profiles of our quasar sample binned in red- 
shift (z < 1.5, and z > 1.5), with median redshift 
z = {0.96,1.96} in the redshift bins, respectively. Since 
our quasar sample is approximately optical flux-limited, 
Malmquist bias causes the higher-redshift quasar sub¬ 
sample to be systematically higher in Mbh and Aboi- 
Due to the nearly redshift-independent nature of the tSZ 
effect and its trends we detect with Mbh and Lboi in 
quasars, the quasar tSZ signal actually increases with 
redshift in our sample. The correlations between red¬ 
shift, Mbhj and Lboi can be isolated by binning our 
quasars in these three parameters simultaneously, but the 
resulting small subsamples causes the signal-to-noise ra¬ 
tio in their stacks to be too low for comparison. We note 
that since the tSZ Compton-y profiles for different quasar 
subsamples are composed of quasars at different angular 
diameter distances and contain contributions from virial- 
ized gas in halos of different masses, the small differences 
in their angular extents are not statistically significant. 


4.3. Implications for Quasar Feedback 

The accretion feedback efficiency of pe ~ 0.05 (or ap¬ 
proximately 5% of the black hole mass) we measure is 
rather discrepant with those needed to reproduce the 
M — a relati on in cosmological ga l axy formation sim - 
ulations (e.g. Di Matteo et al.|[2005 Sijacki et al.||2014 ), 


which generally hnd that ethciencies of r]e ~ 0.005 — 0.015 
are required. We speculate that such large feedback ef¬ 
ficiencies could stem from (1) the incomplete treatment 
of the multi-phase ISM in these simulations and/or (2) 
additional feedback from relativistic quasar jets. 

Regarding (1), since cosmological simulations with 
quasar feedback cannot resolve quasar winds, the sub¬ 
grid models of quasar feedback utilized implicitly assume 
a homogenous ISM with which the quasar wind inter¬ 
acts. However, simulations of the effects of quasar winds 
in more realistic multi-phase ISM environments have 
shown that the hot wind-shocked gas will leak through 
porous low-density regions to galactic scales, allowing the 
higher-density regions t o accrete and produce the ob¬ 


served MBH-g relation (Wagner et a l. 2013 Bourne et 
al. 2014 Zubovas fc Na}^kshinH20Tl). These eftects are 


particularly important in a realistic cosmological context, 
in which filamentary structures of dense infalling gas can 
provide a s ignificant portion of the accreted mass onto 


the SMBH (Costa et al. 2014). Since the quasar feedback 
efficiencies used in simulations are tuned to result in the 
correct observed properties of the simulated galaxies, the 
assumption of a homogeneous ISM can naturally lead to 
the use of small quasar feedback efficiency parameters, 
and correspondingly underpredicted feedback energetics 
in comparison to our observations. 


may be ubiquitous ( 

Fabian et al. 

2014 

) and can con- 

tribute sigmticantly to the teedbac 

c energetics. Most of 


the power in jets - which derives from the angular mo¬ 
mentum of the BH - is eventually converted into thermal 
and kinetic energy, with a maximum power of ~ afMcf 


cretion rate ( 

Blandford & Znajek 

1977 

& McKinney 

2012). Indeed, AcA 

jets 


Tchekhovskoy 
rave been ob- 


timated energetics that reach >10®^erg (Erlund et al. 


2006| Fabian et al.|2014). While it is thought that jets do 
not normally exist during the radiatively efficient quasa r 
state (although see Shaw et ar]|2012 Ruan et aL]|2014[ ), 
they are more comnionly associated with less efficient ac¬ 
cretion states that may have occurred at earlier timesQ 
In order for jets to generate the feedback energetics we 
report, ~10% of the BH mass would have to accreted 
during radiatively inefficient states and the BHs would 
need sizable spin parameters during this previous accre¬ 
tion. Jets would circumvent the feedback constraints 
that derive from reproducing the Mbh ~ cr relation in 
simulations. 

The tentative quasar fe edback detection t hrougj h the 
tSZ effect in WMAP by Chatterjee et al. (20I0| also 
found quasar tSZ signals signihcantly larger (by more 
than an order of magnitude) than expectations from cos¬ 
mological simulations, although no attempt to differenti¬ 
ate between feedback effects and virialized halo gas was 
made in that study. Our detailed consideration of the 
tSZ signal from virialized halo gas in quasars shows that 
it is likely to be subdominant to the additional effects of 
feedback in our observed tSZ signal. Assuming a fiducial 
radiative efficiency oirj = O.I, the large feedback efficien¬ 
cies of e ~ 0.5 of the mean quasar bolometric luminosity 
implied by our results are also consistent with the 0.3- 
I.OLboi commonly observed in the kinetic luminosity of 
near-relativistic quas ar winds as probed through X-ray 


absorption lines (e.g. C hartas et al. 2009 
2009[ Saez et al. 2009 " Chartas et al. 


et al. 


^2014) argued that these large 


Reeves et al. 




T). Chartas 


Bed back efhcien- 
cies are evidence for magnetically-driven quasar winds, 
as opposed to radiative line-driving. Our estimated en¬ 
ergetics are however larger than some other estimates 
of quasa r feedback efficiencies, such a s from FeLoBAL 
quasars ( Faucher-Giguere et al.|[2M^ |. 

The large thermal energy we detect in quasar wind 
bubbles is unlikely to be converted substantially into 
other forms and hence may show up in other tSZ mea¬ 
surements. iScannapieco et al. (2008) considered the 


relic feedback tSZ signal in lower-redshift post-quasar- 
phase galaxies, showing that it is similar in amplitude 
to actively-accreting quasars, and may be detectable by 
stacking Compton-y profiles of low-redshift galaxy pop¬ 
ulations. Differentiating this relic feedback tSZ signal 
from hot virialized gas may be more difficult, since these 
lower-redshift galaxies sit in more massive halos in which 
hot virialized gas contributes more tSZ signal than the 
z ~ 1.5 quasar host halos in our sample. However, semi- 
analytic models place low-redshift IO^M ^t^ black holes in 
IO^^^Mq halos (Wyithe & Loeb 12002 ), for which we 
showed the energetics of virialized halo gas is smaller 
than the feedback energetics, and thus such detections 
of past quasar activity in galaxies remains a possibility. 
Furthermore, at z ~ 0.1 (similar to the galaxy popu¬ 
lations studied in PC13a), the quasar blastwave radius 
of a typical massive galaxy ^10 Gyr after its quasar 
phase is ~5 Mpc, easily resolvable with Planck. 

The Compton-y profile template of extent ~rvir used in 
the matched filter approach of PC13a would not capture 
most of the blastwave Compton-y signal in post-quasar- 
phase galaxies, and this relic tSZ feedback signature may 

^ Since we select quasars that are optically bright and radio 
quiet, it is likely that our systems do not currently have jets. 
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have been missed in these previous investigations. 


5. CONCLUSION 

The effects of quasar feedback through winds have long 
been suspected to shape many of the observed proper¬ 
ties of massive galaxies. We have searched for the effects 
of quasar feedback on ionized gas surrounding quasars 
through the tSZ effect. Using a Planck tSZ Compton- 
y map with explicit corrections for dust contamination, 
we stacked the Compton-y map in regions centered on 
a large sample of SDSS spectroscopic quasars and de¬ 
tected a tSZ signal at >5 (t significance. This stacking 
technique is validated through separate stacks centered 
on SDSS galaxies of different halo masses , from which we 
measured tSZ signals in agreement with Planck Collab- 
oration et al. (2013a). We showed that this tSZ signal 
cannot owe to radio emission or standard parameteriza- 
tions for thermal dust emission, and that it is likely to be 
dominated by feedback effects rather than from virialized 
halo gas. We also find that the signal scales with Mbh 
and Lboi as anticipated for a quasar feedback origin. 

We found the mean angularly integrated Compton-y of 
our sample of quasars to be 3.5 (±0.7) x 10“® Mpc^, and 
estimated the feedback efficiency to be 5% (±1% statis¬ 
tical uncertainty) of the mass of the black hole (although 
with additional systematic uncertainties stemming from 
single-epoch black hole mass estimates. Our measured 
quasar feedback energetics are larger than those used in 
simulations that have been calibrated to reproduce ob¬ 
servations of the Mbh ~ cr relation. If the large quasar 
feedback efficiency we infer is confirmed, we speculate 
that this could be due to (1) the simulations not resolv¬ 
ing the multiphase ISM, which leads to under-predicted 
feedback efficiencies to reproduce the Mbh — cr relation, 
and/or (2) the relic of efficient radio-mode feedback in 
our stack, which taps into the spin energy of the black 
hole. 

Our investigation can be extended using the signih- 
cantly larger samples of quasars and galaxies available in 
current spectroscopic surveys such as the Baryon Oscil- 
lation Spectroscopic Survey (BOSS, Dawson et al.|2013 l 
in SDSS-III, as well as future experiments such as ex- 
tended BOSS (eBOSS) in SDSS-IV a nd the Dark En - 
ergy Spectroscopic Instrument (DESI, Levi et aL||20I3[ ). 
In particular, improved constraints on the quasar halo 
mass distribution from HOD modeling using larger sam¬ 
ples of quasars and galaxies will lead to better account¬ 
ing of the virialized halo gas contribution to the observed 
quasar tSZ signal at low-redshifts, and thus stronger con¬ 
straints on the inferred feedback parameters. Studies of 
correlations between the energetics of quasar feedback 
as probed by the tSZ effect and quasar multi-wavelength 
properties (e.g. dust obscuration, broad absorption lines, 
host galaxy star-formation) can shed light on the driving 
mechanism of the wind and its effects on the host galaxy, 
similar to the trends with Mbh and Lboi that we detect. 
Stacking the tSZ from different lower-redshift galaxies 
populations may also indicate which galaxies have had a 
previous quasar phase. Furthermore, resolved studies of 
the quasar feedback tSZ signal using current higher an¬ 
gular resolution CMB experiments, such as the Atacama 
Cosmolog y Telescope (Hincks et a l. 2010), the South Pole 
Telescope (Carlstro m et al.|201I )7the Ce rro Chajnantor 
Atacama I'elescope (Woody et al. 2012), and Atacama 


Large Millimeter Array (Wootten & Thompson 2009), 
could constrain the spatial extent ot feedback ettects and 
may also be able to better excise troublesome Galactic 
and intrinsic dust emissions. If confirmed, our detection 
of strong quasar feedback will have substantive implica¬ 
tions for our understanding of the role of quasar feed¬ 
back in the evolution of galaxies and of the intergalactic 
medium. 
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